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I  INTRODUCTION 


There  are  several  properties  of  Laplace’s  equation  that  have  important  applications.  One  of  these 
properties  states  that  if  the  magnetic  field  normal  to  the  surface  bounding  a  source-free  region  of  space  is  zero 
eveiywhere  over  that  surface,  then  all  magnetic  field  components  within  the  volume  are  zero  [1],  Another 
related  property  states  that  if  both  orthogonal  components  of  magnetic  field  tangent  to  the  bounding  surface  of  a 
source-free  region  are  zero  everywhere  over  that  surface,  then  all  magnetic  field  components  within  the  volume 
are  once  again  zero  [2],  The  purpose  of  this  paper  is  to  derive  two  analogous  properties  for  special  types  of 
magnetic  field  gradients  called  “directional  derivatives”.  Once  these  properties  have  been  established,  the 
boundary  conditions  ensuring  uniqueness  for  gradient  solutions  inside  the  volume  can  be  determined. 

There  are  many  ways  to  compute  the  gradient  of  a  magnetic  scalar  potential  or  field  component.  In 
general,  the  gradient  of  a  magnetic  scalar  potential,  <p ,  is  expressed  as  Vq>,  and  is  equal  to  the  negative  of  a 

magnetic  field  vector,  H  .  In  Cartesian  coordinates  the  potential  gradient  is  expressed  as 


~d(p  ~d<p  -  d<p 
V®  =7— +  7— +k— 
dx  8y  dz 


(1) 


where  i ,  j,  k  are  unit  vectors  in  the  x ,  y ,  z  directions,  respectively.  In  generalized  curvilinear  coordinates  . 

(1)  must  be  modified  to  account  for  the  curvature  of  the  coordinate  system.  In  an  orthogonal  curvilinear  system 
(OCS),  (1)  is  written  as 


Vp  =  <5, 


dtp 


+  o. 


dm  „  dm 
— — +  77, — — 

h 2&f2 


(2) 


where  <5, ,  a2 ,  a3  are  unit  vectors  in  the  direction  of  generalized  coordinates  ,  £3;  and  where  h] ,  h2 ,  h3 

are  the  metrics  that  scale  the  curvilinear  coordinate  system  under  consideration  [3],  The  rate-of-change  in  each 
component  of  the  vector  in  (2)  is  computed  along  a  curved  coordinate  line  while  holding  the  other  coordinates 
constant. 


To  take  the  process  one  step  further,  the  gradient  of  each  component  in  H  could  also  be  computed.  In 
Cartesian  coordinates,  the  complete  gradient  of  H  is  given  as 


VII  -dH*  -dH*  ldH, 

VHr  =  7 - — +  7 - *-  +  k - 

dx  dy  dz 

(3) 

~dH  .dH  .dH 

VHy=i-f  +  j—^+k-^ 
dx  dy  dz 

(4) 

vu  -dH,  dH  ~  dH 

dx  dy  dz 

(5) 

where  Hx,  Hy,  Hz  are  the  magnetic  field  components  in  Cartesian  coordinates;  while  in  generalized 
coordinates  (reference  3,  p!22,  eq.  7.4)  the  complete  gradient  would  be 


„„  *  dH,  „  dH ,  ,  dH. 

\H.  =77. — 1-77, *-  +  77, - — 

mx  m2  3  m* 


(6) 
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V/72 

VH, 


A  dH2  A  5/7 ,  A  577, 
a. — +o, — +  o, — 

1  vs  2  h2di r2  3 

„  5/7,  A  5/73  A  5/7, 
a , — +«, — +  a, — 

‘Wi  Vl2  3  Via 


(7) 

(8) 


where  77, ,  h2,  H3  are  the  three  components  of  magnetic  field  in  generalized  coordinates.  (It  should  be 
noted  that,  in  general,  the  metrics  A,,  h2,  h3  are  themselves  functions  of  ,  |2 ,  £3 .)  As  was  the  case 

for  the  gradient  of  the  magnetic  potential,  the  gradients  of  the  magnetic  field  vector  components  in  a 
generalized  coordinate  system  give  the  gradient  along  curved  lines  defined  by  the  curvilinear  coordinate. 
This  curvilinear  gradient  is  not  the  boundary  condition  being  considered  in  this  paper. 

In  this  paper,  magnetic  field  gradient  boundary  conditions  will  refer  to  the  directional  derivatives  of 
the  magnetic  field  vectors.  The  directional  derivatives  give  the  rate-of-change  in  magnetic  field  along  a 
straight  line  in  a  specific  direction.  The  directional  derivative  of  the  magnetic  field  vector  H  in  the  direction 
of  a  unit  vector  a  is  given  in  Cartesian  coordinates  by  [4] 

(a  ■  S7)H  =i(a •  V77J  +  j(a  ■  V77v) + k(a  •  V77z)  (9) 


while,  in  generalized  curvilinear  coordinates  the  expressions  are 


[(5- V)//], 


a,  977,  a2  977,  a3  9/7, 
K  9|,  +  h2  9&  +  A3  9£ 


H2  9A,  dh2 

KK  9£2  9#,J 


h3  a \ 
— —  [a,  — - 
Vh  ‘5^ 


dh3 
■a3  — -1 

a I. 


(10) 


[(o- V)77]2 


[(5  •  V)/7]3 


a,  9/72  +  a2  9772  n3  9772 


A,  9£  A2  9|2  A3  9^ 


A 

9A. 


77  3  9A2 

- —  [<3,  — - — n,  — -1 

A,  A,  2  9#3  3  9£2 


*2'*3 

77, 


— 

A,A2  9#, 


3_ 
2 

_  <*}_■. 
'  J 


•2  ^2 
_  a,  9773  n2  9773  9773 

~~h^+Y2l^+T3~d^ 

77,  r  5A3  9A, , 

A,A3 1  3  d^i  9£, 

77 2  r  9A3  9A2 


A2A3  9^2 


9^ 


(11) 


(12) 


where  a, ,  a, ,  a3  are  the  components  of  5  in  curvilinear  coordinates.  The  objective  of  this  paper  is  to 

determine  the  magnetic  field  boundary  conditions  for  directional  derivatives  that  will  guarantee  all  components 
of  the  complete  magnetic  field  gradient  are  uniquely  zero  inside  the  bounded  volume. 
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II  Boundary  Conditions  for  Directional  Derivatives  in  the  Surface-Normal  Direction 

The  standard  derivation  of  the  magnetic  field  properties  of  Laplace’s  equation  will  serve  as  the  bases 
for  proving  the  gradient  cases.  Following  the  procedure  given  in  [1],  applying  the  divergence  theorem  to 
HXVHX  gives 


J  V  •  (HS?Hx)dv  =  §{HXVHX\  hda  (13) 

V  s 

where  da  is  the  differential  area  on  the  closed  surface,  S ,  that  bounds  the  source-free  volume,  V ,  having 
differential  volume  dv  (see  Fig.  1),  and  n  is  the  unit  normal  vector  on  the  bounding  surface  that  points  into  the 
volume.  However, 


V  •  (HXVHX)  =  HXV2HX  +  VHX  ■  VHX 


(14) 


Figure  1.  A  Source  Free  Volume,  V,  bounded  by  the  Surface,  S 


Since  Hx  =  -d<p/dx,  then 
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which  becomes 


<16> 

since  q>  is  a  solution  to  Laplace’s  equation.  (In  general,  the  order  of  differentiation  in  (15)  can  only  be 
reversed  in  Cartesian  coordinates.)  Therefore,  (14)  reduces  to  V-(HXVHX)=  VHx  ■  VHx ,  so  that  (13)  can  be 
written  as 


\f7Hxfdv  =  §Hx(n  ■  VHx)rfa  (17) 

v  s 

Repeating  this  procedure  for  the  y  and  z  components  gives 


and 


\\VHy\dv  =  §Hy{ri-VHy}la 

V  s 


(18) 


j\VHz\2dv  =  §Hz{ri-VHz'yia 

V  s 


(19) 


Adding  (17)  through  (19)  together,  and  rearranging  yields 

/(v^f  +  lv^f+IV^f^v 

V 

=  •  [(«  •  V )FT ]da 


(20) 


The  term  in  brackets  on  the  right  side  of  (20)  is  the  directional  derivative  of  H  in  the  direction  that  is 
normal  to  the  bounding  surface. 

If  all  the  components  of  the  directional  derivative  in  the  direction  normal  to  the  surface  are  zero,  then 
the  right  side  of  (20)  is  zero.  Therefore,  the  integrand  on  the  left  side  is  also  equal  to  zero.  This  means  that  all 
the  magnetic  field  gradient  components  (i.e.  the  complete  gradient  of  H )  are  zero  inside  the  volume. 


It  would  appear  at  first  glance  that  the  zero  requirements  for  all  components  of  the  directional 
derivative  in  the  direction  normal  to  the  surface  applies  only  to  magnetic  field  vectors  expressed  in  Cartesian 
coordinates.  However,  the  result  is  more  general.  If  H  is  expressed  as 


H  =mn+t1ffli+t2Hh  (21) 

where  tx  and  t2  are  unit  vectors  in  the  two  orthogonal  directions  tangent  to  the  surface  (Fig.  1),  then  from  (10) 
through  (12),  the  three  components  of  the  directional  derivatives  in  the  normal  direction  (D  ,D  D  )  are 

v  nn ?  ixn 9  t^n  ^ 

given  by 


D. 


KdS,  KKK, 


(22) 
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D  -  1  dl±  M 
K  S4,  K\ 
z)  _  i  3g<,  n.  dh, 
%  K\  <%, 


(23) 

(24) 


where  hn ,  h,t ,  hh  are  the  metrics  for  the  surface-normal  and  two  orthogonal  tangential  directions, 
respectively.  Therefore,  forcing  Dm  =  Z)  =  Dh„  =  0  ensures  (n  •  V)H  =  0  in  (20). 

When  examining  (20),  it  becomes  clear  that  the  requirement  for  all  three  components  of  the 
directional  derivatives  to  be  zero  may,  in  some  cases,  be  too  rigid.  Expanding  the  integrand  on  the  left  side  of 
(20)  gives 


\VHX\2  +\VHyf  +\VH,\2  =  H2xx  +  H2xy  +  Hl 

+  H2yx+H2yy+H2y!  (25) 

+  HI  +  HI+H2, 

However,  since  H  =H  ,  H  =  H  ,  and  H  =-(H  +  H  ),  it  needs  only  be  shown  that 

|2  +|v/7>|2  =  0  to  prove  that  the  complete  gradient  inside  the  volume  is  zero.  In  fact,  setting  any  two  of 

the  three  terms  on  the  left  side  of  (25)  to  zero  is  sufficient  for  the  complete  gradient  to  be  zero;  which  implies 
that  in  Cartesian  coordinates  only  two  of  the  three  gradient  boundary  conditions  are  required.  There  are  other 
coordinate  systems  in  which  only  two  of  the  three  components  of  the  normal  directional  derivatives  are 
required.  Therefore,  the  general  boundary  condition  specifying  all  three  components  of  the  directional 
derivatives  in  the  surface-normal  direction  is  a  sufficient  condition,  but  not  a  necessary  one. 

Equation  (20)  can  now  be  used  to  prove  the  uniqueness  of  solutions  to  Laplace’s  equations  for  field 
gradient  boundary  conditions.  Assume  there  are  two  possible  solutions  to  Laplace’s  equations  for  field 
gradients  inside  the  volume  V,  Ha  and  Hb .  Let  Hi  =Ha-Hb,  so  that  (h  ■  V)Ha  -  (h  ■  V)Ht  =  0  on  the 
bounding  surface.  (It  is  a  requirement  that  both  solutions  satisfy  the  same  boundary  conditions  on  the 
enclosing  surface.)  Therefore,  since  (n  ■  V)Hd  =  0 ,  |v/f/ 12  +  |VZT*|2  +  |w//|2  =  0  inside  the  volume.  This 
implies  that  the  complete  gradients  of  Ha  and  Hb  are  equal  in  this  region,  proving  uniqueness  of  the  solution. 


Ill  Boundary  Conditions  for  Directional  Derivatives 
in  the  Surface-Tangential  Directions 

The  standard  derivation  of  the  tangential  magnetic  field  boundary  conditions  for  Laplace’s  equation 
will  serve  as  the  bases  for  proving  the  gradient  cases.  Following  the  procedure  given  in  [2],  the  vector  analogue 
of  Green’s  first  identity  is 


J(VxP  -  VxQ  -P  -VxVxQ^Jv 


=  j(P  xV xQ)-hda 

s 


(26) 
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where  as  before  S  is  a  closed  surface  around  the  volume  V .  Let  VH  x  =  VxPx  =  VxQx  (which  implies 

Vx  Vx<2*  =  0),  where  Px  and  Qx  are  given  by  Px  =  a,P*  +  aPx  +  a3Px  and  Qx  =  a,Qx  +  a2Qx2  +  a2Q\ . 

(The  vector  Px  is  shown  to  exist  and  is  defined  in  the  appendix  for  Cartesian  coordinates.)  Rearranging  the 
right  side  of  (26)  gives 


j \VHxfdv  =  jpx-  (VHX  x  h)da 

V  s 

Substituting  n  =  t,  x  i2  in  (27)  and  using  a  vector  identity  from  [4]  gives 


V  S 

Repeating  this  procedure  for  the  y  and  z  components  gives 


J|W,frfv=R-  [ft  • 

v  s 


(27) 


(28) 


(29) 


(30) 


The  terms  inside  the  brackets  on  the  right  side  of  (28)  through  (30)  are  the  directional  derivatives  of 
Hx,  Hy,  Hz  in  the  two  orthogonal  directions  tangent  to  the  bounding  surface.  If  these  six  directional 

derivatives  are  all  zero,  then  the  complete  gradient  within  the  volume  is  also  zero. 

Once  again  it  can  be  shown  that  the  directional  derivative  boundary  conditions  apply  to  curvilinear 
coordinate  systems.  The  six  tangential  directional  derivative  boundary  conditions  from  (28)  through  (30)  can 
be  written  as 


=  VH*)+  •’ k(i,  •  VH, ) 

(f2-V>7 

=  Kh  •  V77x )+  j(t2  •  VHy )+  k{t2  -VHZ) 


Once  again,  if  H  is  expressed  as  in  (21),  then  from  (10)  through  (12),  the  six  components  of  the 
directional  derivatives  in  the  tangential  direction  ( D  ,D  ,D  ,D  ,D  ,D  )  are  given  by 

nh  hh  hh  nh  hh  hh  °  J 


"'1  hn\dt„ 


1  8H,  Ht  dh,  H  dh, 


(33) 

(34) 
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D, 


_  i  mh  Hh  d\ 

\  \K 

_  i  mn  Hh  s\ 
\Hh  KK^n 

__  l  3Hh  Hu  d\ 

8Ht  H  dh.  H.  dh. 


(35) 

(36) 

(37) 

(38) 


where  all  the  terms  have  been  defined  previously.  Therefore,  if  D  =  D  -D  -  D  =  D  -  D  =0,  then 

x  ~  nil  hh  hh  nH  hh  nh 

(irV)H  =(t2-V)H  =  0. 

Using  the  same  arguments  as  for  the  surface-normal  directional  derivatives,  it  can  be  shown  that  all 
six  components  are  not  required  for  a  Cartesian  system,  as  is  the  case  in  some  curvilinear  coordinate  systems. 
Therefore,  the  requirement  that  all  six  tangential  directional  derivatives  be  zero  is  a  sufficient  condition  and  not 
a  necessary  one.  In  addition,  these  six  directional  derivatives  are  sufficient  to  guarantee  a  unique  complete 
gradient  solution  to  Laplace’s  equation. 


IV  CONCLUSION 

Specifying  the  three  magnetic  field  directional  derivatives  in  the  normal  direction  over  a  closed  surface 
around  a  source-free  volume  is  a  sufficient,  but  not  necessary,  condition  to  ensure  a  unique  total  field-gradient 
solution  to  Laplace’s  equation.  Similarly,  specifying  the  six  magnetic  field  directional  derivatives  in  the  two 
orthogonal  tangential  directions  over  a  closed  surface  around  a  source-free  volume  is  also  a  sufficient,  but  not 
necessary,  condition  to  ensure  a  unique  total  field-gradient  solution  to  Laplace’s  equation.  Although  the  proofs 
in  this  paper  were  conducted  specifically  for  magnetic  fields,  the  directional  derivative  boundary  conditions 
should  apply  to  any  field-gradient  boundary  value  problem  in  physics.  In  addition,  it  may  be  possible  to  use 
this  theory  to  develop  finite  element  techniques  that  use  directional  derivatives  as  input  to  boundary  value 
problems  and  compute  unique  solutions  for  the  complete  gradient  inside  a  volume. 


7 


NSWCCD-TR- 1999/007 


REFERENCES 


Paris,  D.  T.  and  F.  K.  Hurd,  Basic  Electromagnetic  Theory,  New  York:  McGraw-Hill  Inc  1969 
p.165. 


Stratton,  J.  A.,  Electromagnetic  Theory,  New  York:  McGraw-Hill,  1941,  p.256. 

Tai,  C.  T,  Generalized  Vector  and  Dyadic  Analysis,  2nd  Edition,  New  York  and  Oxford:  IEEE  Press 
and  Oxford  University  Press,  1997,  p.  62. 

Morse,  P.  M.  and  H.  Feshbach,  Methods  of  Theoretical  Physics,  Parti,  New  York'  McGraw-Hill 
1953,  p.  33. 


8 


NSWCCD-TR- 1999/007 


APPENDIX 

In  this  appendix  it  will  be  shown  that  a  vector  Px  can  be  found  such  that  VHX  =Vx?'.  The 

magnetic  field  H  can  be  represented  as  the  curl  of  a  vector  potential  A  (H  =  V  x  A  )  and  can  be  written  in 
Cartesian  coordinates  as 


dy  dz  dz  dx  dx  dy 


where  Ax,  Ay,  A.  are  the  components  of  A  .  From  (Al),  the  expression  for  VHX  is 


VHx  =  i(- 


d2A, 


dydx  dzdx 


f,d2A 

+  k( 


d2Ay 

dzdy 

d2A 


) 


dydz  dz‘ 


-) 


Next,  VxP1  is  given  by 


VxPx=i  ( 


dP[ 

dy 


.-5£)+aE 

dz  ’  dz 

-  dPx 
+  k( — — 


dx 


dPx 


dx  dy 


) 


Equating  like  components  between  (A2)  and  (A3)  gives  the  three  equations 

dPx  dPx.  _  d2Az  d2 Ay 

dy  dz  dydx  dzdx 

dP*  dPx  _  d2A2  d2Ay 

dz  dx  dy2  dzdy 

dP>  dPx  _  d2Az  d2  Ay 
dx  dy  dydz  dz 2 

From  (A4)  it  is  clear  that  Px  and  Px  can  be  represented  as 


(Al) 


(A2) 


(A3) 


(A4) 

(A5) 

(A6) 


(A7) 

(A8) 


since  the  order  of  differentiation  can  be  changed  in  Cartesian  coordinates.  What  remains  is  to  show  that  the 
selection  of  Px  and  Px  in  (A7)  and  (A8)  are  consistent  with  (A5)  and  (A6).  Placing  (A8)  into  (A5)  and  taking 

the  partial  derivative  with  respect  to  y  gives 
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a2p;  _  d  ,d2Az  |  d2Az)  d2Ay 
dydz  dy  dx 2  dy 2  dy2dz 


(A9) 


However,  Hx  is  also  a  solution  to  Laplace’s  equation  as  shown  by  (15)  and  (16)  of  the  main  text,  and 
can  be  written  from  (Al)  as 


ck2  dy  dz  dy2  dy  8z  * 


+  31(54.  _fi)  =  0 

dz  dy  dz ' 


which  can  be  rearranged  to  give 


a  ,a24  d2A,  d2A: . 

dyKdx2  dy2  dz2’ 

d  ,d2A  d2A  d2A 

—  — ( - —  H - — H - £-) 

dz  dx 2  dy2  dz2 


Using  (All)  with  (A9)  and  reducing  produces 

gg  _  a  gj, 

Sy&  a?  av2  <3z2  ayaz 


(A  10) 


(All) 


(A12) 


Placing  (A7)  into  (A6)  and  taking  the  partial  derivative  with  respect  to  z  gives  an  expression  that  is 
identical  to  (A12).  Therefore,  the  choice  of  p*  and  P*  given  in  (A7)  and  (A8)  are  consistent  with  (A4) 
through  (A6). 


In  this  appendix,  it  has  been  shown  that  P *  exists  such  that  VHX  =Vxf'.  In  addition,  the 

expressions  for  the  three  Cartesian  components  of  P*  have  been  determined.  Similarly,  it  can  be 
shown  that  P y  and  Pz  exist  such  that  VH  =V  xPy  and  VH,  =  VxP. 
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